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ABSTRACT 

Radio-wave scintillation observations reveal a nearly Kolmogorov spectrum of density fluctuations in 
the ionized interstellar medium. Although this density spectrum is suggestive of turbulence, no theory 
relevant to its interpretation exists. We calculate the density spectrum in turbulent magnetized plasmas 
by extending the theory of incompressible magnetohydrodynamic (MHD) turbulence given by Goldreich 
& Sridhar (1995) to include the effects of compressibility and particle transport. Our most important 
results are as follows. 

(1) Density fluctuations are due to the slow mode and the entropy mode. Both modes are passively 
mixed by the cascade of shear Alfven waves. Since the shear Alfven waves have a Kolmogorov spectrum, 
so do the density fluctuations. 

(2) Observed density fluctuation amplitudes constrain the nature of MHD turbulence in the interstellar 
medium. Slow mode density fluctuations are suppressed when the magnetic pressure is less than the gas 
pressure. Entropy mode density fluctuations are suppressed by cooling when the cascade timescale is 
longer than the cooling timescale. These constraints imply either that the magnetic and gas pressures 
are comparable, or that the outer scale of the turbulence is very small. 

(3) A high degree of ionization is required for the cascade to survive damping by neutrals and thereby to 
extend to small lengthscales. Regions that are insufficiently ionized produce density fluctuations only on 
lengthscales larger than the neutral damping scale. These regions may account for the excess of power 
that is found on large scales. 

(4) Provided that the thermal pressure exceeds the magnetic pressure, both the entropy mode and the 
slow mode are damped on lengthscales below that at which protons can diffuse across an eddy during 
the eddy's turnover time. Consequently, eddies whose extents along the magnetic field are smaller than 
the proton coUisional mean free path do not contribute to the density spectrum. However, in MHD 
turbulence eddies are highly elongated along the magnetic field. From an observational perspective, the 
relevant lengthscale is that transverse to the magnetic field. Thus the cut-off lengthscale for density 
fluctuations is significantly smaller than the proton mean free path. 

(5) The Alfven mode is critically damped at the transverse lengthscale of the proton gyroradius, and 
thus cascades to smaller lengthscales than either the slow mode or the entropy mode. 

Subject headings: MHD — turbulence — ISM: kinematics and dynamics 



1. INTRODUCTION 

Diffractive scintillations of small angular-diameter radio 
sources indicate that the interstellar electron density spec- 
trum on lengthscales 10^ — 10^° cm is nearly Kolmogorov; 
i.e. r.m.s. density fluctuations across a lengthscale A are 
nearly proportional to A^/^. They also establish that there 
are large variations in the amplitude of the density spec- 
trum along different lines of sight. 

Rickett (1977, 1990) and Armstrong, Rickett, & Span- 
gler (1995) review the observations of diffractive scintilla- 
tion and their interpretation. They also discuss refractive 
scintillations and dispersion measure fluctuations, which 
probe density fluctuations on scales larger than the diffrac- 
tive scales. Non-diffractive measurements tend to indicate 
that the Kolmogorov spectrum extends to much larger 
scales. However, we focus primarily on diffractive mea- 
surements because they are much more sensitive. 

Density fluctuations that obey the Kolmogorov scal- 
ing occur in homogeneous subsonic hydrodynamic turbu- 



lence. They are due to the entropy mode, a zero-frequency 
isobaric mode whose density fluctuations are offset by 
temperature fluctuations. Since subsonic turbulence is 
nearly incompressible, the velocity fluctuations follow Kol- 
mogorov's scaling. To a good approximation, the entropy 
mode is passively mixed by the velocity field, so it also con- 
forms to the Kolmogorov spectrum.^ Density fluctuations 
in the Earth's atmosphere, which cause stars to twinkle, 
obey the Kolmogorov scaling. They arise from the passive 
mixing of the entropy mode. 

The electron density spectrum in the interstellar 
medium cannot be explained by hydrodynamic turbu- 
lence.^ Because the medium is ionized, magnetic effects 
must be accounted for. This is evident since the length- 
scales probed by diffractive scintillations are smaller than 
the collisional mean free paths of both electrons and pro- 
tons. If the magnetic field were negligible, freely streaming 
plasma would wipe out density fluctuations at diffractive 
scales. In the presence of a magnetic field, electrons and 



^ Density fluctuations due to the Reynolds stress scale as A^/^ . 

^ Charge neutrality is maintained on diffractive scales, so electron density fluctuations include compensating fluctuations in the density of 
positive ions. 
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protons are tied to fieldlines at the scale of their gyroradii. 
For typical interstellar field strengths, these gyroradii are 
smaller than the diffractive scales. A magnetic field thus 
impedes the plasma from streaming across fieldlines, and 
allows the turbulent cascade and the associated density 
fluctuations to reach very small scales across the fieldlines 
before dissipating. Therefore, a theory for compressible 
turbulence in magnetized plasmas is required to explain 
the observed density spectra. Our objective is to develop 
this theory. 

Until now, the only description of density fluctuations 
in interstellar plasmas was by Higdon (1984, 1986). These 
papers, while prescient, preceded a theory for incompress- 
ible MHD turbulence, and therefore did not account for 
the full dynamics of the cascade. We compare Higdon's 
theory with ours in §10. 

Our compressible theory extends the theory of incom- 
pressible MHD turbulence given by Goldreich & Sridhar 
(1995) by including a slightly compressible slow mode and 
a passive entropy mode. We also consider kinetic effects: 
on sufficiently short lengthscales, the mean free paths of 
the particles are significant, and the equations of compress- 
ible MHD must be modified. This is especially important 
for damping. 

In a future paper II, we apply the theory developed here 
to estimate amplitudes for density fluctuations produced 
in supernova shocks, HII regions, stellar winds, and ther- 
mally unstable regions. These are then compared to scat- 
tering measures observed along different lines of sight. 

Before considering compressible turbulence, we discuss 
incompressible MHD turbulence, focusing on issues that 
are important for the compressible case. 

2. INCOMPRESSIBLE MHD TURBULENCE 

Goldreich & Sridhar (1995, 1997) propose a picture of 
the dynamics of incompressible strong MHD turbulence 
and describe the power spectra of Alfven waves, slow 
waves, and passive scalars. We extend their picture to 
cover additional features such as the parallel cascades of 
both slow waves and passive scalars. Throughout this pa- 
per, "parallel" and "transverse" refer to the orientation 
relative to the "local mean magnetic field", which is the 
magnetic field averaged over the scale of interest. Our 
discussion of incompressible MHD turbulence, while some- 
what lengthy, is important for understanding the extension 
to compressible turbulence that follows. 

Consider a uniform unperturbed plasma with an cm- 
bedded magnetic field. Turbulence is excited at the MHD 
outer scale, Lmhd, by random and statistically isotropic 
forcing, with r.m.s. velocity fluctuations and r.m.s. mag- 
netic field fluctuations (in velocity units) which are com- 
parable to the Alfven speed, va-^ 

As the turbulence cascades from the MHD outer scale to 
smaller scales, power concentrates in modes with increas- 
ingly transverse wave vectors. The inertial range velocity 
spectrum applies to lengthscales below Lmhd but above 
the dissipation scale. It is anisotropic and is characterized 
by 

= VA I J 1 , (1) 

^ The forcing fluctuations may also be less va, in which case Lmhd 

to VA- 
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(2) 



the inertial range magnetic field spectrum is identical. 
Here is the lengthscale transverse to the local mean 
magnetic field, v\j^ is the r.m.s. velocity fluctuation across 
A^, and A|| is the lengthscale parallel to the local mean 
magnetic field across which the velocity fluctuation is vxj^ ■ 
We interpret Ay as the elongation along the magnetic fleld 
of an "eddy" which has a size A_l transverse to the mag- 
netic field; it is not an independent variable, but is a 
function of A^. Deep within the inertial range, where 
A^ <C Lmhd, eddies arc highly elongated along the mag- 
netic field: A|| X±. In the following subsections, we 
explain the physics underlying the spectrum, and consider 
some of the implications. 

2.1. Alfven Wave Spectrum 

Arbitrary disturbances can be decomposed into Alfven 
waves and slow waves. The Appendix summarizes the 
properties of these waves in the more general case of com- 
pressible MHD. In incompressible MHD, Alfven waves and 
slow waves are usually referred to as shear- Alfven waves 
and pseudo- Alfven waves, but the former designation is 
more convenient for making the connection with compress- 
ible MHD. 

Our understanding of the MHD turbulence is based on 
two facts: (i) MHD wave-packets propagate at the Alfven 
speed either parallel or antiparallel to the local mean mag- 
netic field; and (ii) nonlinear interactions are restricted to 
collisions between oppositely directed wave-packets. These 
facts imply that in encounters between oppositely directed 
wave-packets, each wave-packet is distorted as it follows 
field lines perturbed by its collision partner. A wave- 
packet cascades when the fieldlines that it is propagating 
along have spread by a distance comparable to its trans- 
verse size. 

Alfven waves have quasi two-dimensional velocity and 
magnetic field fluctuations which are confined to planes 
perpendicular to the local mean magnetic field. As their 
more complete name shear- Alfven implies, they dominate 
the shear of the mapping of planes transverse to the local 
mean magnetic field produced by field line wander. Thus, 
Alfven waves control the dynamics of MHD cascades; slow 
waves may be ignored when considering the dynamics of 
Alfven waves. 

In strong MHD turbulence the cascade time of an Alfven 
wave-packet is comparable to its travel time across the 
parallel length of a single oppositely directed Alfven wave- 
packet of similar size. Goldreich & Sridhar (1995) refer to 
this balance of timescales as "critical balance" . It relates 
the parallel size of a wave-packet, A|| , to its transverse size, 
A^. Wave-packets of transverse size A_l cascade when the 
fieldlines they follow wander relative to each other by a 
transverse distance A_l. Critical balance implies that this 
occurs over a parallel distance A||. 

The Alfven wave spectrum is given by equations (1) 
and (2), with referring to the velocity fluctuations 
of the Alfven waves. It is deduced from two scaling ar- 
guments: (i) Kolmogorov's argument that the cascade 
time tx^ ~ X±/vx^^ leads to an energy cascade rate, 

would be defined as the lengthscale at which the fluctuations extrapolate 
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v\^/tx^ — ■'^Aj^/'^-Li which is independent of lengthscale; 
and (ii) the "critical balance" assertion that the linear 
wave period which characterizes the Alfven waves in a 
wave-packet is comparable to the nonlinear cascade time 
of that wave-packet, i.e. va/ ^\\ — ^/t\^- 

Before considering slow waves in MHD turbulence, we 
discuss two topics that are governed by the dynamics of 
Alfven waves only: eddies and passive scalars. 

2.2. Eddies 

Because of their transverse polarization, Alfven waves 
are responsible for the wandering of magnetic fieldlines. 
A snapshot of wandering fieldlines is shown in Figure Al. 
Each of these fieldlines passes through a localized region 
of size A_L in one plane transverse to the mean magnetic 
field. Away from this plane the bundle of fieldlines di- 
verges due to the differential wandering of the individual 
lines. At a second plane, the bundle's cross sectional area 
has approximately doubled. Critical balance implies that 
the distance to this second plane is comparable to the par- 
allel wavelength which characterizes the bimdlc, A||. As 
the bundle spreads, other fieldlines, not depicted, enter 
from its sides. In general, the neighboring fieldlines of 
any individual field line within a region of transverse size 
A_L change substantially over a parallel distance of order 
A||. It is natural to think of A|| as the parallel size of an 
"eddy" that has transverse size A_l. Two eddies with the 
same transverse lengthscale that are separated by a par- 
allel distance greater than their Ap incorporate different 
fieldlines, and hence are statistically independent. Eddies 
are distinct from wave-packets. The former arc rooted in 
the fiuid whereas the latter propagate up and down mag- 
netic fieldlines at the Alfven speed. 

Aside from their anisotropy, eddies in MHD turbulence 
are similar to those in hydrodynamic turbulence. They 
are spatially localized structures with characteristic veloc- 
ity fluctuations and lifetimes. The r.m.s. velocity dif- 
ference between two points is determined by the small- 
est eddy that contains both. Different eddies of a given 
size are statistically independent. The three-dimensional 
spectrum for r.m.s. velocity fluctuations across transverse 
lengthscales Aj_ and parallel lengthscalcs A|| is 

V. . - X / (^^/Lmhd)1/3 , for A|| « A|| 

There is negligible additional power within an eddy on 
parallel lengthscales smaller than Ay, so for Ay <C A||(A_l), 
^'Aj_,A|[ = ^'Aj_- For A|| » A||(A^), the smallest eddy that 
contains both A_l and Ay has a transverse lengthscale A^ 
which satisfies A||(A^) = A||. The velocity fluctuation of 
this eddy is obtained by solving this equation for X'^ (cq. 
[2]), and inserting this A'^ in equation (1). Contours of 
the three-dimensional spectrum are plotted in Figure A2. 
Each contour represents eddies of a characteristic size. 

Maron & Goldreich (2001) give the three-dimensional 
spectrum in Fourier-space. Since eddies that are separated 
by more than Ay arc statistically independent, the power 
spectrum at a fixed transverse wavenumber k±_ is indepen- 
dent of the parallel wavenumber Aiy in the corresponding 

region of Fourier-space, i.e. where A;|j~^ > A||(fcJ^). 

The turbulent cascade is generally viewed as "proceed- 
ing" from larger eddies to smaller eddies as this is the di- 



rection of energy transfer. However, smaller eddies cascade 
many times in the time that a large eddy cascades. This is 
particularly important in turbulent mixing. Consider the 
evolution of two fluid elements whose initial separation is 
larger than the dissipation scale. On cascade timcscale 
, their transverse separation will random walk a dis- 
tance A_L as the result of the cascade of smaller eddies. 
Therefore, on a timescale comparable to an eddy's cas- 
cade time, the transverse locations of its component fluid 
elements — whose sizes may be considered to be compara- 
ble to the dissipation scale — are completely randomized. 
Moreover, since mixing at the dissipation scale causes 
neighbouring fluid elements to be rapidly homogenized, 
transverse smoothing of the eddy occurs on the timescale 
that it cascades. Rapid transverse mixing in MHD tur- 
bulence is similar to the more familiar isotropic mixing in 
hydrodynamic turbulence. 

2.3. Passive Scalar Spectrum 

A passive scalar, cr, satisfies the continuity equation, 
(d/dt + V ■ V)(J = 0, and docs not affect the fluid's evo- 
lution. It could represent, for example, the concentration 
of a contaminant. We consider the spectrum of a passive 
scalar mixed by the Alfven wave cascade. These consid- 
erations are important for our subsequent investigation of 
compressible turbulence. They arc also helpful for under- 
standing the slow wave spectrum. We discuss the passive 
scalar spectrum both in the inertial range and also below 
the scale at which the Alfven wave spectrum is cutoff. 

2.3.1. Passive Scalar Spectrum in the Inertial Range 

As we show in this subsection, the transverse spectrum 
of the passive scalar in the inertial range is 

cx X'f , (4) 

where A_l is the lengthscale transverse to the local mean 
magnetic field and dx^ is the r.m.s. fluctuation in the 
passive scalar across A^. The passive scalar parallel spec- 
trum is the same as the Alfven wave parallel spectrum 
given in equation (2), where Ay is now to be interpreted 
as the lengthscale parallel to the local mean magnetic field 
across which the passive scalar fluctuation is (7x± ■ 

Mixing of the passive scalar is due to Alfven waves. Slow 
wave mixing is negligible. This is because transverse veloc- 
ity gradients are much larger than parallel ones in MHD 
turbulence. Thus Alfven waves, whose velocity fluctua- 
tions are perpendicular to the magnetic fleld, are much 
more efii'ective at mixing than slow waves, whose velocity 
fluctuations are nearly parallel to the magnetic field. The 
transverse cascade arises from the shuffling of fieldlines as 
Alfven waves propagate through the fluid. 

The transverse spectrum (eq. [4]) follows from the 
Kolmogorov-like hypothesis that the cascade rate of the 
"energy" in the scalar field is independent of lengthscale; 
i.e. crx^/tx± is constant, where tx^ is the passive scalar 
cascade time, which is assumed to be proportional to the 
cascade time of Alfven waves. Comparing this with the 
constancy of the kinetic energy cascade rate, v'^^/tx^, we 
conclude that ax^ oc vx^, which implies equation (4). A 
similar argument holds for the cascade of a passive scalar 
in hydrodynamic turbulence (e.g., Tennekes & Lumley 
1972). 
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The parallel cascade of a passive scalar is more subtle. 
It niight appear that a passive scalar cannot cascade along 
fieldlines since, neglecting dissipation, both the scalar and 
the magnetic field are frozen to the fluid, and thus the 
scalar must be frozen to fieldlines. In that case there cer- 
tainly could not be a parallel cascade. If the scalar were in- 
jected on large scales, then fluctuations with smaller wave- 
lengths along the magnetic field would not be generated. 
However, dissipation cannot be neglected. It is an essential 
part of MHD turbulence, as it is of hydrodynamic turbu- 
lence. For example, the description of turbulent mixing in 
§2.2 depends crucially upon small scale dissipation. 

Perhaps the best way to understand the parallel cas- 
cade is to consider mixing on the transverse lengthscale 
Ax within two pianos which arc perpendicular to the local 
mean magnetic field, and which are separated by a parallel 
distance larger than Ay. Velocity fluctuations within the 
two planes are statistically independent. This is evident 
because a bundle of fieldlines cannot be localized within a 
transverse distance A_l over a parallel separation greater 
than A||. Even a pair of fluid elements, one in each plane, 
that are initially on the same fieldline are mixed into two 
regions with different values of passive scalar concentra- 
tion. It follows that the parallel cascade of the passive 
scalar also obeys equation (2). 

2.3.2. Passive Scalar Spectrum Below Alfven Wave 

Cut- Off 

A passive scalar cascade may extend below the trans- 
verse scale at which the MHD cascade is cut-off. Mixing 
on these scales is driven by fluid motions at Acutoff which 
results in a scale-independent mixing time equal to the 
cascade time at Acutoff- This yields 

(Ja^ = constant , Aj_ < Acutoff • (5) 
A similar argument applies in hydrodynamic turbulence. 
Tennekes & Lumley (1972) call this regime in hydrody- 
namic turbulence the "viscous-convective subrange" . 

2.4. Slow Wave Spectrum 

The slow wave spectrum is the same as that of the 
Alfven waves. It is given by equations (1) and (2), with 

V\j^ referring to the velocity fluctuations of the slow waves. 
This is a consequence of the similar kinematics of slow 
waves and shear Alfven waves and the fact that both are 
cascaded by shear Alfven waves. 

Slow waves obey the same linear wave equation as 
Alfven waves, and to lowest nonlinear order they travel up 
and down the local mean magnetic field linos at the Alfven 
speed just as Alfven waves do. However, the dynamics of 
the MHD cascade is controlled by the Alfven waves (Gol- 
dreich & Sridhar 1997; Maron & Goldreich 2001) because 
their velocity and magnetic field fluctuations arc perpen- 
dicular to the local mean magnetic fleld, whereas those of 
the slow waves are nearly parallel to it. Since perpendic- 
ular gradients arc much larger than parallel ones in the 
MHD cascade, Alfven waves arc much more effective at 
mixing than are slow waves. Hence, Alfven waves cascade 
both themselves and slow waves, whereas slow waves cas- 
cade neither.* 

The transverse mixing of the slow waves by Alfven waves 
is analogous to the mixing of a passive scalar. As discussed 

* We are assuming that Alfven and slow waves have comparable str< 



in §2.3, a passive scalar assumes the same inertial range 
spectrum as that of the velocity field which is responsible 
for its mixing. Thus, equation (1) is also applicable to the 
velocity fluctuations of the slow waves. 

Similarly, the parallel cascade of slow waves is analogous 
to the parallel cascade of a passive scalar. Since Alfven 
waves cascade in the time they move a distance Ay, slow 
waves separated by this distance are independently mixed. 
Thus, equation (2) also applies to slow waves. There is, 
however, a conceptual difference between the parallel cas- 
cades of the passive scalar and of the slow mode. In the 
absence of dissipation a passive scalar is frozen to field- 
lines, whereas slow mode wave-packets travel along them 
at the Alfven speed. A passive scalar has a parallel cascade 
bcc;ausc Alfvcnic fluctuations are statistically independent 
within two transverse planes frozen in the fluid and sep- 
arated by A||. The parallel cascade of slow waves occurs 
because two transverse planes separated by A[| and trav- 
elling at the Alfven speed in the same direction experience 
uncorrelated sequences of distortions suffered as a restilt of 
interactions with oppositely directed Alfven waves. Nev- 
ertheless, these two requirements are both satisfied in the 
MHD cascade when the transverse planes are separated by 
a distance greater than Ay, and so the passive scalar and 
the slow mode have the same parallel spectrum. Whereas 
passive scalar mixing is due to eddies, slow wave mixing 
may be thought of as due to "travelling eddies" , i.e. eddies 
which travel up and down the magnetic field at the Alfven 
speed. 



2.5. Numerical Simulations 

Numerical simulations offer some support for the above 
description of incompressible MHD turbulence. Those 
by Cho & Vishniac (2000b) support both equations (1) 
and (2), and those by Miiller & Biskamp (2000) sup- 
port equation (1). However, although the simulations of 
Maron & Goldreich (2001) support equation (2), they yield 

v\j_ oc A^'^^ instead of equation (1). Because the sim- 
ulations of Maron & Goldreich (2001) arc stirred highly 
anisotropically, whereas those of Cho & Vishniac (2000b) 
and Miiller & Biskamp (2000) are stirred isotropically, it is 
not clear whether these disparate results conflict. Maron & 
Goldreich (2001) speculate that the discrepancy between 
their spectrum and the scaling prediction of Goldreich & 
Sridhar (1995) results from intermittoncy. In any case, we 
expect that the physical picture of a critically balanced 
cascade, which underlies Goldreich & Sridhar's description 
of MHD turbulence, remains valid. Even if the spectrum 

is proportional to A^''^, we expect that the results of this 
investigation which assumes a spectrum proportional to 

1/3 

A^ — would not be significantly altered. 

The simulations of Maron & Goldreich (2001) confirm 
that a passive scalar has the same transverse spectrum as 
that of Alfven waves, although both are proportional to 
A^^^. They also indicate that the parallel cascade of a 
passive scalar conforms to equation (2). 

Maron & Goldreich (2001) present results from simula- 
tions of the interaction between oppositely directed slow 
and Alfven waves; the slow waves cascade whereas the 

IS at a given lengthscale. 
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Alfvcn waves do not. They also compute the spectrum 
of slow waves in a simulation of MHD turbulence and find 
that both its transverse and longitudinal behavior matches 
that of the Alfven wave spectrum. 

2.5.1. Kolmogorov Constants 

Scaling arguments do not yield values for the "Kol- 
mogorov constants" , the order-unity multiplicative con- 
stants of the spectrum. However, they can be obtained 
from simulations. We define them such that equation (1) 
remains valid, i.e. we take Lmhd to be the separation at 
which the r.m.s. velocity difference is equal to, or extrapo- 
lates to, va- Two Kolmogorov constants, M|| and Mt, are 
needed in this paper: 

M|| ^ ^ All = M,| UI/^Lmhd^/^ (6) 



In these definitions, Ax and Ay are the inverses of the 
wavcnumbcrs transverse and parallel to the local mean 
magnetic field, and i^j. is the cascade time of waves with 
transverse wavenumber AJ^. From the numerical simula- 
tions of Maron & Goldreich (2001), 

M|| ~ 3.4 and ~ 1.4 . (8) 
Because these simulations yield a transverse spectrum 
which is proportional to A^^''^ instead of A^^'"* , the resulting 
"Kolmogorov constants" are not truly constant. 

3. COMPRESSIBLE TURBULENCE: OVERVIEW 

Our primary concern is interstellar scintillation, which 
is affected by electron density fluctuations on very small 
scales, typically 10^ — 10^° cm for diffractivc scintillation. 
In the remainder of this paper we calculate the spectrum of 
density fluctuations which results from compressible tur- 
bulence in magnetized plasmas with parameters appropri- 
ate to the interstellar medium. Throughout, we consider 
plasmas that have more ions than neutrals, and that have 
1 < /3 < oo, where /3 is the ratio of the thermal pressure 
to the mean magnetic pressure: /3 = Snp/B^ = 2c^/v\. 
Here, ct = {p/pY^"^ is the isothermal sound speed, va 
is the Alfven speed, p is the thermal pressure, B is the 
magnetic fleld strength, and p is the mass density.^ The 
incompressible limit corresponds to /3 = oo. 

On the lengthscales that we consider, compressible 
MHD is a good approximation for the dynamics of the 
ionized interstellar medium. Kinetic effects, where impor- 
tant, may be accounted for by simple modifications to the 
MHD equations. Therefore, we turn our attention to tur- 
bulence in compressible MHD. 

The turbulent velocity spectrum in compressible MHD 
is approximately the same as the turbulent velocity spec- 
trum in incompressible MHD, because the Alfven mode 
remains incompressible in a compressible medium, and the 
slow mode is only slightly compressible. Thus, the velocity 
spectrum for both of these modes is given by equations (1) 
and (2). The Appendix summarizes the properties of the 
relevant modes in compressible MHD. 

There are two additional modes in compressible MHD 
which are not present in incompressible MHD. One of these 

^ We briefly discuss plasmas with /3 < 1 in §11. 



is the fast mode. However, as long as /? > a few, the fast 
mode is essentially a sound wave. Its phase speed is larger 
than the phase speed of either the Alfven mode or the 
slow mode, and so the fast mode does not couple to them. 
Thus, we ignore the fast mode. An analogous approxima- 
tion is often made in subsonic hydrodynamic turbulence, 
where sound waves have little influence on the turbulent 
cascade. 

The second mode that is present only in compressible 
MHD is the entropy mode. This is a zero-frequency mode 
with unperturbed pressure, and with density perturbation 
offset by temperature perturbation (see Appendix). The 
entropy mode exists when the turbulent motions are adia- 
batic; the entropy mode is then the only source of entropy 
fluctuations. In hydrodynamics, simple scaling arguments 
show that the entropy mode does not affect the dynamics 
of the fluid — i.e. the fluid obeys the incompressible equa- 
tions of motion — provided that fluid motions are subsonic 
and that density fluctuations are smaller than the mean 
density; these arguments carry over directly to MHD (e.g. 
Higdon 1986, and references therein). Furthermore, since 
the entropy of any fluid element is conserved in the inertial 
range of adiabatic turbulence, the mixing of the entropy 
mode is identical to the mixing of a passive scalar (see eq. 
[4]), yielding the transverse spectrum 

sx^ oc , (9) 

where saj^ is the r.m.s. entropy fluctuation across A_l . The 
parallel spectrum is given by equation (2). Analogously, in 
hydrodynamic turbulence the effect of the entropy mode 
on turbulent motions can often be neglected, and the en- 
tropy mode is mixed as a passive scalar (Monin & Yaglom 
1971). 

Based on the above discussion, there are two sources of 
density fluctuations on small scales: the slow mode and the 
entropy mode. Since the slow mode density perturbation 
is proportional to its velocity perturbation, and since the 
entropy mode density perturbation is proportional to its 
entropy perturbation, both yield a Kolmogorov spectrum 
of density perturbations: 

nx^ (X XY^. (10) 

In the remainder of this paper, we investigate the den- 
sity spectrum in more detail. We calculate the density 
spectrum on diffractive lengthscales for given values of the 
number density, n, and outer scale. Lout- Typically, in 
regions of the interstellar medium which are relevant for 
scintillation, 1 cm^'^ < n < 100 cm^^. Values of Lout 
are more uncertain, though it is likely that Lout is typi- 
cally within a few orders of magnitude of 1 pc. It is also 
plausible that the value of Lout is related to that of n. In 
paper II, we discuss in much more detail the values of these 
parameters in turbulent regions of the ionized interstellar 
medium, e.g. in supernova shocks, HII regions, and stellar 
winds. 

A naive guess for the resulting density fluctuation at the 
lengthscale A^ is = n(A_L/Lout)^^'^, which we refer to 
as the flducial spectrum. However, there are a number 
of physical effects which suppress the small-scale density 
spectrum in the interstellar medium relative to the fiducial 
spectrum. 
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kinetic energy density is sufficiently small that it is com- 
parable to the mean magnetic energy density: ?^Lmhd ~ '^A' 
which implies that 



JMHD 



Lout/3 



-3/2 



(11) 



Below this scale, the turbulent kinetic energy density is 
smaller than the magnetic energy density, and the theory 
of compressible MHD turbulence is applicable. Note that 
the effects of large-scale velocity fields can be neglected 
when considering fluctuations on smaller scales, because a 
large-scale velocity field can be eliminated by a Galilean 
transformation. However, since large-scale magnetic fields 
cannot be transformed away, they affect the dynamics of 
small-scale eddies. 

Hydrodynamic turbulent motions on scales slightly 
larger than Lmhd have speeds ^ va- Thus they couple 
to, and efficiently excite, Alfven waves and slow waves 
whose phase velocities u/k k VA^kg/k) arc of comparable 
magnitude. Provided (3 > 1 and M < I, fast waves, which 
have u) /k <^ Cs > VA, are weakly excited at Lmhd- In what 
follows, we neglect fast waves. 

Alfven waves excited at Lmhd cascade to smaller scales, 
and these cascading Alfven waves in turn cascade slow 
waves. The spectra of both of these cascades are given by 
equations (1) and (2), so that 



A I \i/3 



(12) 



When considering scintillation observations, it is the 
transverse not parallel — lengthscale which is relevant. 
Each line of sight crosses regions with different orientations 
of the local mean magnetic field. The observational effects 
of the parallel lengthscale arc washed out if the orientation 
varies by an angle greater than a tiny ratio: A_l/A|| at the 
diffractive scale. Since the variation in angle due to eddies 

1/3 

on scales larger than the diffrative scale is A^/A|| cx A^^' , 
it increases with scale, so these eddies render the parallel 
lengthscale unobservable. As a result, wo frequently refer 
to the transverse lengthscale as, simply, the lengthscale. 

Because of the multitude of special lengthscales in- 
volved, we organize the discussion by decreasing length- 
scale. We begin at the outer scale, and proceed down 
to the smallest scales relevant to interstellar scintillation, 
while considering the effects of the Alfven mode, the slow 
mode, and the entropy mode simultaneously. Since most 
of the relevant lengthscales in the interstellar medium have 
a similar dependence on the background density — they 
decrease with increasing density — the ordering of length- 
scales is fairly universal. The lengthscales which we con- 
sider are the outer scale, the MHD scale, the cooling scale, 
the coUisionless scale of the neutrals, the collisionless scale 
of the ions, and the proton gyro scale. We conclude with 
a summary of the most important effects. Table A. 3 lists 
the lengthscales that are used in this paper. 

4. THE OUTER SCALE AND THE MHD SCALE 

As a model for the excitation of the turbulence, we con- 
sider plasma which is stirred on an outer scale Lout with 
velocity fiuctuations that are of order the sound speed: 
^Lout ^ Cg. The generalization from this case of Mach 
1 turbulence to subsonic turbulence with arbitrary Mach 
number M. < I is trivial: Lout would be interpreted as 
an effective outer scale at which velocity fiuctuations ex- 
trapolate to Cg. However, we focus on the Mach 1 case 
because it is probably the most relevant for interstellar 
scintillation. 

If, initially, the strength of the magnetic field is negli- 
gible, then random fieldline stretching amplifies the mean 
field (Cho & Vishniac 2000a). It is uncertain both how 
quickly the magnetic field is amplified, and whether its en- 
ergy density is amplified until it approaches equipartition 
with the turbulent kinetic energy density. If it does reach 
equipartition within a few outer-scale crossing times, then 
/? ~ 1 would be appropriate for Mach 1 turbulence. Recall 
that /3 is the ratio of thermal pressure to magnetic pres- 
sure, or equivalently, (3 = 2c^/v\. We generally assume 
that the mean magnetic field can approach equipartition 
with the gas pressure, so we think of /3 as close to but a 
little larger than unity, although we leave its exact value 
unspecified. 

Provided (3 is larger than unity and M. Ion scale Lout, 
the turbulent kinetic energy dominates the mean mag- 
netic energy on scales just below Lout- Thus the cascade 
is hydrodynamic on these lengthscales, and the velocity 
fiuctuations are given by Kolmogorov's isotropic scaling: 

^ Cs(A_L/Lout)^/^- 

The kinetic energy density decreases towards smaller 
scales. At a critical scale — which we denote Lmhd — the 

® For definiteness, we consider plasma that is thermally stable; deviations of the temperature from its equilibrium value decay in a cooling 
time, tcool- 



for both modes, even on scales smaller than Lmhd- 

We have not yet discussed density fluctuations on the 

scales which have been considered in this section. It is 
more convenient to do so in the following section. 

5. THE COOLING SCALE 

On the lengthscales in the interstellar medium that we 
have considered thus far, the radiative cooling time is 
shorter than the eddy turnover time.^ Consequently, in 
photoionized regions, turbulence on these lengthscales is 
expected to be isothermal (Sridhar 1998; Goldreich 1998; 
Higdon & Conley 1998). This has only a marginal effect on 
the turbulent dynamics described in the previous section, 
because isothermal Alfven waves are identical to adiabatic 
Alfven waves, and isothermal slow waves are only slightly 
different from adiabatic slow waves. However, there is no 
entropy mode in isothermal turbulence. As a consequence, 
small-scale density fluctuations may be significantly sup- 
pressed. There arc two possible solutions to this "cool- 
ing catastrophe": either (i) the outer scale is extremely 
small, small enough that the turbulence at the outer scale 
is nearly adiabatic; or (ii) there arc significant density fluc- 
tuations associated with the slow mode. However, in the 
latter case, the mean magnetic field must be amplified al- 
most to equipartition with the gas pressure, so that /? ~ 1. 
Either of these two solutions would place stringent con- 
straints on the nature of the turbulence which is responsi- 
ble for observed density fluctuations. In the following, we 
consider the cooling scale in more detail. 

There is a critical scale for the turbulence, which we call 
the "cooling scale", Lcooi- Above this scale the turbulence 
is isothermal, and below it the turbulence is adiabatic. We 
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assume that Lcooi < Lmhd throughout this paper, except 
in §5.2 where we consider the case in which this inequality 
is reversed. The cooling scale is where the eddy turnover 
time, A_l/waj^, is equal to the cooling time, tcooi; i-e., with 

VX^ ~ Cs(A_L/Lout)^/^, 



Lmit >' 



-L'out ^ J^out ' 

As we will discuss in detail in paper II, most plausi- 
ble astrophysical sources of interstellar scintillation have 
Cs^cooi ^ Lout, and thus Lcooi ^ Lout- In general, the 
kinetic power per unit mass which is dissipated in Mach 
1 turbulence is ~ c^/th^ut ~ c^/Lout, where ^Lout is the 
eddy turnover time at the outer scale. The thermal power 
per unit mass which is required to keep the gas hot is 
^ Cg/tcooi- Hence, to avoid the conclusion that the cool- 
ing scale is smaller than the outer scale, one must require 
that more energy go into turbulent motions than into heat. 
As we will sec in paper II, this requirement is difficult to 
satisfy. Nonetheless, in isothermal shocks, and in plasmas 
which are thermally unstable, the two powers are com- 
parable, and the cooling scale is comparable to the outer 
scale. 

As an example, we consider an HII region. Photoion- 
ized plasma is thermally stable. Heating is primarily due 
to photoionizing photons, and cooling is primarily due to 
electron impact excitation of metal line transitions (e.g., 
Spitzcr 1978).'' A characteristic temperature for plasma 
photoionized by hot stars is T ~ 8,000K, which implies 
that the speed of sound is Cg ~ 10 km/sec. The cool- 
ing time is tcooi ^ 20, 000( cm~^/n) years, where n is the 
number density of electrons. Therefore, 



0.2 



(^) 



pc 



(14) 



A plausible value for Lout might be the radius an HII re- 
gion (the "Stromgren radius"), which is ~ 70( cm~^/n)^/^ 
pc (Spitzcr 1978). Therefore, with typical values of n in 
HII regions (1 cm~'^ < n < 100 cm~^), Lout is significantly 
larger than Cstcod, and so the cooling scale is significantly 
smaller than the outer scale. 

The fact that the turbulence is isothermal on large scales 
has important implications. Had cooling been ignored, i.e. 
had it been implicitly assumed that the cooling time is infi- 
nite, then one would have calculated the density spectrum 
as follows: there should be density fluctuations of order 
unity on the outer scale, implying excitation of entropy 
modes on the outer scale which arc passively mixed by the 
Alfvenic turbulence to small scales. This would yield the 



fiducial density spectrum, nx^ ~ n(A^/Lo 



U/3 



However, if the turbulence is isothermal on large scales, 
the above calculation is incorrect: the entropy mode does 
not exist under isothermal conditions. In the following, 
we consider separately two cases: first the low-/3 case (/? < 
Lout/csicooi), and then the high-/? case (/3 > Lout/csicooi)- 
In each case, wc show that the small-scale density fluc- 
tuations are substantially smaller than those predicted by 



the fiducial spectrum unless a relatively extreme condition 
holds: either /3 1 or Lout ~ Cgtcooi- 

5.1. Density Fluctuations Below the Cooling Scale: 

cooi 

With regards to interstellar scintillation, the crucial dif- 
ference between high-/? and low-/? turbulence lies in the 

compressibility of the slow mode. Since the sum of ther- 
mal pressure and magnetic pressure vanishes for the slow 
mode, the mode's density perturbation satisfies nx^/n 

Vx^/p - ~ [3-h^x^/vA ^ p-^l^vx^/c, (see 

Appendix). So, assuming that the energy in slow waves is 
comparable to that in Alfven waves, 

nx^ 1 / A_L y/3 

n v^VLout^ 



A± <L 



MHD 



(15) 



from the slow mode, both above and below the cooling 
scale. Thus, slow waves can produce density fluctuations 
which are not much smaller than the fiducial spectrum if 
(3 is not much larger than unity. Comparing with equa- 
tion (25) below, the contribution of the slow mode to the 
density fluctuations at Lcooi exceeds the contribution from 
Reynolds stresses provided /? < Lout/csicooi- 

The compressibility of the slow mode also yields density 
fluctuations associated with the entropy mode. The rea- 
son is as follows. Entropy fluctuations are associated with 
isothermal compressible waves. Those due to isothermal 
slow modes are passively mixed by Alfven waves. Mixing 
due to Alfven waves with wavelengths smaller than the 
cooling scale takes place in less than the cooling time. It 
produces a spectrum of entropy modes for A^ < Lcooi that 
gives rise to a density spectrum similar to the one given in 
equation (15). Nonetheless, there is some damping of en- 
tropy fluctuations when they are mixed from scales larger 
than the cooling scale to scales smaller than the cooling 
scale. We expect that the damping is of order unity. In this 
paper, however, we do not quantify the amount of damp- 
ing more precisely, because it is a difficult calculation; it is 
more difficult than the damping of the slow mode, which 
we calculate in the following section. 

5.1.1. Slow Mode Damping at the Cooling Scale 

Although there is negligible damping of the slow mode 
on lengthscales which are cither much larger or much 
smaller than the cooling scale, there is some damping on 
lengthscales which are comparable to the cooling scale. 
As we show in this section, the damping of the slow wave 
cascade is small, because the slow mode is not very com- 
pressible. In the limit that /3 ^ oo, the slow mode is 
incompressible, and the damping disappears. In the fol- 
lowing, we calculate the damping to first order in 1 //?. 

First, we calculate the damping rate of a slow wave of 
a fixed wavenumber. From the Appendix, the dispersion 
relation of the slow mode, in the limit that fcii/fc <C 1, is 



VAk\\ 



(16) 



In this case, there is another lengthscale, the "photoionization scale", which we do not consider because its effects are unimportant for 
scintillations. This is the scale at which the eddy turnover time is comparable to the recombination time. It is slightly larger than the cooling 
scale because the recombination time is approximately five times larger than the cooling time (Spitzer 1978). On scales larger than the pho- 
toionization scale, the turbulence is in photoionization equilibrium, whereas on scales smaller than this, the ionization fraction of a fluid element 
is conserved. Nonetheless, the temperature of the gas is only weakly dependent on ionization fraction: metal-line cooling by free electrons is 
exponentially dependent on temperature, and thus acts as a thermostat which is only slightly affected by the ion or neutral density. Thus, the 
photoionization lengthscale does not play an important role in the density spectrum. 
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The notation which we use in this section is the same as 
the notation which we define in tlic Appendix, except that 
here we label the axis parallel to the mean magnetic field 
by II instead of by z. Note that we define = p'/p'> 
where primes denote Eulcrian perturbations. The lin- 
earized equation of energy conservation is 



iujTs' = 



V 



cool 



(17) 



where s' is the perturbed entropy per particle. Combin- 
ing this equation with the following monatomic ideal gas 
relations 



yields 



p p T 3 3p' 



~2 _ 2 2 - i5ujtcoo\ 

p' 2 - iSwicooi 



(18) 



(19) 



where cf. = p/p = (iv\/2 is the square of the isother- 
mal sound speed. On large scales, wtcooi 1, and c is 
the isothermal sound speed; on small scales, ujtcooX S> 1, 
and c is the adiabatic sound speed. The damping rate is 
obtained by substituting the above c (eq. [19]) into the 
dispersion relation, which gives, to lowest order in 1//3, 



-VAknt. 



I ''COOl 



VAk\\ /3[1 + (5wAfc||icooi/2)^ 



(20) 



whore uii is the imaginary part of w. Only on Icngthscales 
whore i'a^ii^cooI ~ 1 is this ratio non-nogligiblc. 

To obtain the total decrement in the amplitude of the 
slow mode through the cooling scale, we solve a kinetic 
equation obtained by balancing the slow mode fc-space en- 
ergy flux with the loss-rate of slow mode A:-space energy 
density due to damping; i.e. 



dk 



k 



(21) 



^ n,^ is the transverse wavenumber, t\_^ is the 



where = k 
cascade time of a slow wave with this wavenumber, and 
is the velocity perturbation of the slow mode; since 
fc_L ~ k, wc drop the ± subscript on k. 
kinetic equation as follows: 



We rewrite the 



where Mt = VAk\\t\^ ~ 1.4 is a Kolmogorov constant (see 
§2.5.1). Substituting equation (20) into equation (22) and 
integrating, using k\\ oc /c^/^ (eq. [2]) and t^j. oc k~'^/^, 
yields the net damping through the cooling scale: 



5.2. Density Fluctuations Below the Cooling Scale: 

P > Lout/Csicool > 1 

In this subsection only, wc assume that /3 > Lout/cs^cooi, 
which implies Lmhd < Lcooi- In this case, the isothermal 
hydrodynamic turbulent cascade extends from the outer 
scale down to the cooling scale. Because the entropy mode 
does not exist under isothermal conditions and hydrody- 
namic tm'bulcncc is incompressible to order v\^/ca, to 
this order there are no density fluctuations on lengthscales 
larger than Lcooi- However, pressure fluctuations due to 



Reynolds stresses, 
sity and entropy fluctuations, 

n p \ Cs 



create second order den- 



^cool 



(24) 

Those entropy fluctuations couple to the entropy mode at 
the cooling scale. On scales below the cooling scale, the 
entropy mode is mixed as a passive scalar yielding 



"A, 



A, \V3 



/ J-icool \ 

\ Lniit / 



Lcooi 
2/3 



^ Lcooi ' 



Lout 

Cstcool\l/2/ A_L 



1/3 



1/3 



1/311 

J_ J|Ai<Leool 



= exp 



cxp 



(_Mt f°° _ 
\ 13 Jo 1 

) 



VAk\\tcool 



dk^ 



1 + (5uAfc||tcooi/2)2 k 



10/3 



~ 1 - 1.3//3 



(23) 



to lowest order in 1/(3. Thus for (3 slightly larger than 
unity, slow mode damping can be ignored. 



J , , T I , A_L < Lcooi • (25) 

Hydrodynamic turbulence mixes the entropy mode down 
to Lmhd- At this scale, the hydrodynamic motions turn 
into Alfven and slow waves, so the entropy mode contin- 
ues to be mixed below Lmhd, and equation (25) remains 
valid. Comparing with equation (15), the density fluc- 
tuations associated with the entropy mode exceed those 
associated with the slow mode when /3 > Lout/csicooi- 

The entropy mode's density spectrum (eq. [25]) is 
smaller than the flducial spectrum by the small factor 
(cs^cooi/Lout)"'^^^- As will be seen in paper H, it yields 
density fluctuations that are too small to explain much 
of the observed interstellar scintillation for most plausible 
values of Lout, such as those given by Stromgren radii of 
HII regions. There fire two ciltcriiatives; cither Lout 

is not 

much larger than Cgtcooi: or /3 < Lout/cstcooi- 

6. THE COLLISIONLESS SCALE OF THE NEUTRALS 

In the remainder of this paper, we consider lengthscales 
which are comparable to, or smaller than, the mean free 
paths of the various particles; thus, kinetic effects must be 
considered. 

In this section, wo calculate the damping of the tur- 
bulence by neutrals. We assume that the neutral den- 
sity is smaller than the electron density, i.e. tiat < n, 
and that most of the ions are protons. The neutral 
particles which we consider are hydrogen atoms and he- 
lium atoms. The relevant lengthscales are the hydrogen 
and helium mean free paths for collisions with protons: 
Lh = 5xlQ^^{cm-^/n) cmandLne = 1.5x ^^^(cm-^/n) 
cm (Banks 1966). The hydrogen mean free path is signif- 
icantly smaller than the helium mean free path because 
its collisions with protons are due to resonant charge ex- 
change. 

There are two main results for this section. 

• (i) Only if the neutral fraction is sufficiently small 
(eq. [53] , below) can the cascade survive on length- 
scales smaller than Lh- If the neutral fraction does 
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not satisfy this condition, then ah three modes — 
Alfven, slow, and entropy — are damped at the same 
lengthscale. 

• (ii) Regions where the cascade does not survive con- 
tribute to an excess of density fluctuations on large 
scales. This might explain observations which de- 
tect an excess of power in large-scale density fluctu- 
ations. 

The organization of the calculation is as follows; first, 
the change in wave; frequency due to neutrals is calculated. 
Second, we consider the effect of the frequency change on 
the cascade. Third, we consider the cases in which the 
neutrals are hydrogen atoms and in which they are helium 
atoms. Although helium atoms have a larger mean free 
path, we consider them after neutral hydrogen atoms be- 
cause they are of lesser importance. Finally, we consider 
regions where neutral damping terminates the cascade. 

A similar, though less detailed calculation is performed 
by Goldreich & Sridhar (1995) for the case in which the 
neutrals are hydrogen atoms. 

6.1. Frequency Change 

Consider an incompressible wave, either an Alfven wave 
or a slow wave in the incompressible limit. Although the 
slow wave is, in fact, slightly compressible, this does not 
have a large effect on the final result. We calculate the 
frequency change of an incompressible wave with fixed k: 
/cy along the mean magnetic field and transverse to the 
mean magnetic field, where ^ fc_L ~ fc. 

We define v to be the mean velocity of the protons av- 
eraged over the proton distribution function. Thus v sat- 
isfies the equation of motion for the Alfven wave or for 
the slow wave derived in the Appendix. We then calculate 
the force that these protons exert on the neutrals. Since 
this force is the same as the force of the neutrals on the 
protons, inserting it into the equation of motion for the 
protons yields the change in frequency. 

We denote the perturbed distribution function for the 
neutrals, i.e. perturbed to first-order in vn, after being 
Fourier-transformed in space and time, by SfN{k,u!,VN), 
where Vat is the velocity of the neutrals. The evolution of 
the neutrals is determined by the Boltzmann equation: 



-iuj S/n + i{k ■ VN)SfN = C 



(26) 



where C is the Fourier-transformed, perturbed, collision 
integral of the neutrals with the protons. Neutral- neutral 
collisions are negligible relative to neutral-proton colli- 
sions. The collision integral is simplified by making the 
approximation that, on the timescale that a neutral col- 
lides with protons, S/n is driven towards a Maxwellian 
with mean velocity v: 



{v ■ Vn) riM cxp{-mNvjf/2T) 



5f. 



N 



, (27) 



where the first term in the square brackets is generated by 
expanding a Maxwellian distribution function with mean 
velocity v to linear order in v and retaining only the per- 
turbation, and vn,p is the frequency with which a given 
neutral atom collides with protons. With the above ap- 
proximation to the collision integral, the solution to the 



Boltzmann equation reads 

r., {v ■ Vn) UN exp{-mNV%/2T) 

OjN = 



l'N,p 



T/niN (27rT/mjv)^/^ vn p — ico + ik ■ vn ' 

(28) 

Next we verify that the perturbed neutral number den- 
sity vanishes. 



/ 



SfNd^VN = 



UnWn^p 



(27r)3/2(T/mAr)5/2 



dvN\\VN\\ exp(-mArw^||/2T) 



/ 



d'v 



exp{-mNV%^/2T) 



(29) 

'^N,p — lOJ + ik ■ Vn±- 

= , (30) 

since the first square-bracketed integral vanishes by anti- 
symmetry. For the purpose of evaluating the above in- 
tegrals, we use _L and || for the perpendicular and par- 
allel projections oi vn relative to v, and not relative 
to the magnetic field as in the rest of this paper; i.e. 
vn\\ = {v ■ vn)v/v^ and vn± = vn — vn\\- We also 
use the incompressibility relation k ■ v = 0, which implies 



that k ■ v 



N\\ 



0. 



With the neutral distribution function given by equa- 
tion (28), the force per unit volume of the protons on the 
neutrals may be calculated as 



(31) 



F = TUN J VN{-iuj6fN + i{k ■ VN)5fN)d^VN 
_ mNTlNVUN^p 



(27r)3/2(T/miv)5/2 
j dvN\\v%\\ex.p{-mNV%\\/2T) 

{—iui + ik ■ vn±) exp(— mivt''|^^/2rj 



/ 



d vn±- 



vn p — + ik ■ vn± 



^2) 



mNnNVVN,p 
(27rr/mjv)i/2 

{—icu + iks) exp(— mivs2/2T) 



ds- 



mNnNVVN,p 



^N,p — -|- iks 



7ri/2ni 

=° It 



l + it 



exp [- (t/Hi -hn2/ni)2]dt 



(33) 



(34) 



The second equality follows after replacing the overall mul- 
tiplicative Vn (= vn\\ + vnl) by i'Ar||, because the vni. 
term vanishes by antisymmetry of the integral over Vn\\- 
In the second equality, the first bracketed integral is equal 
to {2'k)^/'^{T /mN)'^^'^ ■ The second bracketed integral is a 
double integral; the integral over the component of t;jv_L 
which is perpendicular to k is equal to {2ttT /ttin)^^'^ ■ The 
remaining integral is over the component of vn_l which is 
parallel to fc; we label this component s in the third equal- 
ity. Finally, the fourth equality follows from the change of 
variables t = {ks — u))/ i'N,p, and from the following defini- 
tions 



Hi = 



k / 2T \ 1/2 



Ho 



vn,p ^ ruN I 

uj VAk\\ 



\mN^ 



vn,p vn,p 



(35) 



(36) 
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where nip is the proton mass. The dimensionlcss parame- 
ter Hi is the number of wavelengths a neutral crosses be- 
fore colliding, and the dimensionless parameter 112 is the 
number of waveperiods a neutral travels before colliding. 
Since k ^ fc|| and (3 > 1, Hi ^ II2 in the inertial range of 
the MHD cascade. 

Performing the integral for the real and imaginary parts 
of F to lowest order in 112/111 yields 

F = mNnNVVN,p{^ - 21^^51(111) , (37) 



where 



ff(ni) = 1 - Viexp(nr^)erfc(nr')/ni 

in? , for Hi < 1 
1 , for Hi » 1 • 



(38) 
(39) 



Since — F is equal to the force of the neutrals on the pro- 
tons, we insert —F into the equation of momentum conser- 
vation for the protons. More precisely, we add —iF/ {nipn) 
to the right-hand side of equation (A2) in the Appendix. 
We label the resulting change of frequency caused by the 
presence of neutrals as Au. Assuming that Au -C VAk\\, 
it follows that 
Aw 



VAk 



mr,n 



1 

nf 



2U2 



.9(ni) 



(40) 



6.1.1. 



Discussion of the Frequency Change 

The physical interpretation of this frequency change is 
straightforward. Recall that k ^ k\\, so that the wave- 
length, 27r/fc, is nearly identical to the wavelength trans- 
verse to the magnetic field, 2TT/k±. It is convenient to 
define the neutral mean free path, Ljv, as follows: 



Hi _ ct rzmp \ 
k I'N.p ^ rnN ^ 



Ct /2m„\i/2 



(41) 



where ct — [p/ p)^^'^ is the isothermal sound speed. 

The imaginary part of Aw is the negative of the damping 
rate, and is given by 

Wi = -- VN,pg{kLN) (42) 



1 rriNnN 



for fcLjv <C 1 
for fcLjv 1 



(43) 



The above expression may be explained as follows (Gol- 
dreich & Sridhar 1995). For wavelengths larger than the 

neutral mean free path, fcLjv < 1, the neutrals are nearly 
locked to the protons. They damp the MHD waves by 
transferring the protons' momentum across a neutral mean 
free path. The damping rate increases with decreasing 
wavelength provided fcL^v < 1- For wavelengths smaller 
than the neutral mean free path, the neutral atoms are 
effectively freely streaming. This leads to a damping rate 
which is independent of wavelength. 

The real part of the frequency change, Awr, is 



AUr 
VAk\\ 



1 

(khN) 



for /cLjv -C 1 
for fcLjv '> 1 



(44) 



Its physical significance is apparent. For small wavelengths 
the neutral atoms are effectively freely streaming, the mo- 
tion of the protons nearly decouples from the motion of 
the neutrals, and there is negligible real frequency change 
associated with the presence of the neutrals. For large 
wavelengths the neutrals are locked to the protons, so the 
mass density of the fluid which participates in the waves 
is larger than that of the protons. Since the Alfven speed, 
and hence the wave frequency, is inversely proportional to 
the square root of the mass density, i.e. uj o: va p~^^'^, 
this results in a decrease of the frequency by a fractional 
amount which is equal to one-half of the mass ratio of 
neutrals to protons. 

6.2. Effect of Neutrals on the Turbulent Cascade 
6.2.1. Effect ofwi on the Cascade 
The turbulent cascade is quenched if the damping 



(X 



rate exceeds the eddy cascade rate or, equivalently, 
> 1. From equation (43) and k 

is seen that for small wavenumbers 

and for large wavenumbers |a;i|/w^fc|| cx 



if 

A:2/3, it 
\/vAki\ cx k^/^, 
^-2/3_ Therefore 



|a;,|/?;^fc|| is a maximum for transverse wavelengths, k~^, 
comparable to the neutral mean free path, and decreases 
for both larger and smaller wavelengths. The requirement 
that the cascade survive damping by neutrals is then that 
|w,|/t;Afc|| < 1 for k~^ comparable to the neutral mean free 
path, or equivalently, nNi'N,p/n < VAk\\ at this scale (eq. 
[43]). Since nNiy^^p/n is the rate at which a given proton 
collides with neutrals, the cascade survives to small scales 
if the wave frequency, and hence the cascade rate, at the 
scale of the neutral mean free path is faster than the rate 
at which a proton collides with neutrals. 

To obtain the total decrement in the amplitude of both 
the Alfven mode and the slow mode through the damping 
scale, we solve the kinetic equation for the cascade. * As in 
§5.1.1, the kinetic equation is obtained by balancing the 
fc-space energy flux with the loss-rate of fc-space energy 
density due to damping: 



— ^ = 2w,: 



V 



(45) 



dk 

X — fc is the transverse wavenumber, and 
t\^ is the cascade time. The kinetic equation can be 
rewritten as follows: 



where Xj^ = k 
the 



d 



V 



1 ^-L 

— — - m — — 
dink txj^ 



2Mt 



VAku 



(46) 



where Mt is a Kolmogorov constant (see §2.5.1). Before 
integrating this equation, we re-express the damping fre- 
quency in terms of the relevant lengthscales. Lout and Ln- 



VAk\\ 



1 mNriN 



2M|| 
1 



VAk\\ 

(fcLjv) 



(47) 



\/2M|| 



V rUp J n 



k 

/ Lout \ 



(48) 



g(fcLjv) 

(fcL;v)2/3 



(.49) 



* The calculation of the amplitude decrement in this section is similar to that given in §5.1.1 for slow mode damping at the cooling scale. The 
principal difference, aside from the fact that the imaginary part of the frequency is different, is that here we must account for the decrease 
in the cascade rate, txj_ , caused by damping. This is not necessary when treating slow mode damping at the cooling scale because there the 
Alfven modes which control the cascade rate are undamped. 
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The first equality above follows from equation (42). The 
second equality follows from the Kolmogorov constant 
M|[ = VAk\\/vx^k (see §2.5.1). The third equality fol- 
lows from the velocity spectrum, vx_i^ = CT(fcLout)~^^'^ (eq. 
[12]), and from the definition of the neutral mean free path 
(eq. [41]). For convenience, we define Lout in this section 
such that = ct when = Lout- 

Integrating the kinetic equation (cq. [46]), and using 
^Aj. ~ Aj./i'Aj^, yields the net decrement due to damping 
by collisions with neutrals: 

[^Ai/Ay^]|A^»LN 



= exp 



: exp 



2Mt 



dk 







VAk\\ k 

1/2 



V TOr) / n V Lat / J 



Lat 



after substituting from equation (49). The numerical pref- 

actor in the second equality, 0.2, follows after inserting the 
values of the Kolmogorov constants (eq. [8]) and the value 
of the integral 



g{x)x ^/^di 



1.1 



(51) 



which was integrated numerically; see equation (38) for 
the definition of g{x). 

For the cascade to continue to small scales, the right- 
hand side of equation (50) cannot be very small. This 
can be viewed as an upper limit on the neutral fraction: 
jiN jn < (Ljv/Lout)^^^- Recall that the condition for the 
cascade to reach small scales is that the cascade rate at 
Ltv be faster than the rate at which a proton collides with 
neutrals. If the neutral fraction is too large, then so is the 
proton collision rate, and the cascade is quenched. More- 
over, for a fixed Ljv, a large value for Lout implies that the 
cascade time at Ljv is large, and hence that the cascade is 
more susceptible to damping by neutrals. 

We postpone consideration of the damping due to hy- 
drogen and helium atoms until after we evaluate the effect 
of on the cascade. 

6.2.2. Effect of AoJr on the Cascade 

We can picture the cascade as proceeding from large 
scales to small scales. As it crosses the scale of the neu- 
tral mean free path (cq. [44]), the effective Alfven speed, 
and hence the real part of wave frequency, increases by the 
fraction ruNriN /2mpn. Consequently, the cascade time de- 
creases by the same amount. To the extent that the flux 
of energy in Alfven waves from large scales to small scales 
is constant in the turbulent cascade, v\^/t\^^ = constant, 
where txj^ is the cascade time. Therefore, a fractional de- 
crease of TUNnM /'i.mpn in the cascade time causes a decre- 
ment in the small-scale Alfvcnic velocity perturbation by 
the fractional amount mNUN / ^mpTi. This decrement is 
in addition to that due to damping. Moreover, it applies 
to the slow mode and the entropy mode as well as to the 
Alfven mode. Because riN/n must be small for the cas- 
cade to pass through the scale of neutral damping, this 
decrement is also small, and we ignore it from here on. 

6.3. Neutral Hydrogen Atoms 

When the neutrals are hydrogen atoms we set njv = nu, 
i'N,p = i^H.p, Ljv = Lh, and mjv = trip. Collisions between 



hydrogen atoms and protons are due to resonant charge ex- 
change. From equation (41), with the value of y-a_,p taken 
from Banks (1966), 

Lh = 5x 10^3/^^) cm (52) 

at 8 X 10^ K; the temperature dependence of Lh is very 
weak. 

For the cascade to survive damping by neutral hydrogen 

atoms, the right-hand side of equation (50) cannot be very 
small. This sets an upper limit on the neutral fraction of 

nH< / exlO^^cm y/3 

n ~ VLout^* V(n/cm-3)Lout>' ' ^ ' 

The real part of cj/w^fc|| is larger on scales below Lh than 
it is on scales above Lh by the fractional value nii/2n. 

6.4. Neutral Helium Atoms 

Although helium has a lower abundance than hydrogen, 
it has a higher ionization potential. Therefore, in some 
regions helium might comprise the majority of neutrals. 
When the neutrals are helium atoms we set njv = nue, 

I'N.p = I'Hc.p, Ljv = Lhc, and ~ Amp. From equation 
(41), with the value of i^He,p taken from Banks (1966), 

-3 s 



^He 



1.5 X 10^^ 



(^) 



cm 



(54) 



at 8 X 10-'' K, and Lho oc T^/^. 

To place an upper bound on damping by helium atoms, 
we assume that most of the helium is neutral, i.e. riHe ~ 
O.It?. Then the decrement is given by the right-hand side 

of equation (50): 

r /\i/3i 



^1/31 



Only if 

Lout ^ 2 X 10^^ (cm ^/n) cm could the cascade be 
terminated at the scale of the helium mean free path. If 
there are not many neutral helium atoms, or if the outer 
scale is not sufficiently large, then damping at Lhc may be 
neglected, and the cascade extends at least to the scale of 
the hydrogen mean free path. 

Collisions of neutral helium atoms with singly ionized 
helium ions might also be significant. Although there are 
fewer helium atoms than protons, He''-He"'" collisions have 
a larger cross-section than He°-proton collisions because 
they are due to resonant charge exchange. Nonetheless, 
given the cosmic abimdance of helium, the mean free path 
of neutral helium due to Hc"-He+ collisions cannot be sig- 
nificantly smaller than that due to He"-proton collisions 
regardless of the ionization fraction of helium. 

The real part of w/uaAh increases by 2n]ie/n below the 
scale Lhb, or 0.2 if most of the helium is neutral. 

6.5. // Neutrals Damp the Cascade 

In subsequent sections, we consider regions in the inter- 
stellar medium where there are too few neutrals to damp 
the cascade. However, there are almost certainly many re- 
gions where the neutrals do damp the cascade. We discuss 
these regions here. 

Suppose both Alfven waves and slow waves are damped 
by neutrals. What happens to the entropy waves? If un- 
damped, they would be mixed down to smaller scales by 
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Alfven waves at the neutral damping scale. The resulting 
density spectrum would be = constant (see eq. [5]). 

1/3 

Because it is flatter than nx^ oc , regions in which 
damping by neutrals truncated the Alfven cascade might 
be important contributors to small scale density fluctu- 
ations. However, the fact that they would contribute a 
density spectrum different from that which is observed 
suggests that the entropy wave cascade is not more re- 
sistant than the Alfven wave cascade to damping by neu- 
trals. Indeed that is the case. Recall that the condition 
for the truncation of Alfven and slow wave cascades is that 
each proton collide with at least one neutral atom during 
one wave period at the neutral damping scale. Under this 
condition, the neutrals would damp the entropy waves by 
conducting heat across them. 

Although regions in which the cascade is damped by 
neutrals do not contribute small scale density fluctuations, 
they may still be significant. Observational cvidc;iic;c; indi- 
cates that there is more power in density fluctuations on 
large scales, 10^^ — 10^^ cm, than would be predicted by 
extrapolating from small scales, 10^ — 10"'^'' cm, with the 
Kolmogorov scaling. See, for example, Lambert & Rick- 
ett (2000) for a review of the observations. Perhaps this 
excess arises in regions where the cascade is damped by 
neutrals. 

We complete this section by briefly considering and then 
rejecting the possibility that a turbulent cascade truncated 
at the neutral damping scale might be regenerated on a 
much smaller scale due to stirring by eddies at the damp- 
ing scale. Although the ratio of the damping rate to the 
wave frequency decreases below the scale of the neutral 
mean free path, the absolute damping rate approaches a 
constant value. Provided the cascade is truncated by neu- 
tral damping, this rate is larger than the stirring rate and 
the cascade cannot be regenerated. 

7. THE COLLISIONLESS SCALE OF THE IONS 

If the cascade survives the neutral collisionless scales, 
then, proceeding to smaller scales, the next scale of im- 
portance is the ion collisionless scale. This scale is set by 
the mean free path of protons to collide with other protons: 

Lmfp = 6 X 10"(^^) cm (56) 

at a temperature of 8,000K (Braginskii 1965). Since the 

electron and proton gyroradii (see eq. [66] for protons) 
are very small compared to the proton mean free path, 
the electrons and protons are tied to magnetic fieldlines. 
Therefore, when considering the collisionless effects of elec- 
trons and protons, the relevant lengthscales are those par- 
allel to the magnetic field. Since turbulent eddies are 
highly elongated along the magnetic field, their transverse 
lengthscales are much smaller than their parallel ones (eq. 
[2]). It is the transverse lengthscale which is relevant when 
considering observations of the density spectrum, because 
each line of sight averages over regions with different mag- 
netic field orientations. In this section, we show that both 
the slow mode and the entropy mode are cut off at the 
lengthscale where the parallel eddy size is comparable to 
the proton mean free path. The lengthscale where the 
density spectrum is observed to cut off, i.e. the trans- 
verse lengthscale, is therefore significantly smaller than the 
mean free path. 



Before discussing the damping of the slow mode and the 
entropy mode, we consider two larger lengthscales, which 
are set by the effects of the electrons. Although the be- 
haviour of the cascade at these lengthscales is interesting, 
it is shown to be unimportant for our purposes. 

Throughout this section, we neglect numerical factors of 
order unity, such as the factors which are associated with 
kinetic corrections to the fluid equations (given in Bragin- 
skii 1965) and the Kolmogorov constants. However, we 
retain the dependences on f3, which is assumed to satisfy 
/? > 1, and on the ratio of the proton mass to the electron 
mass: nip /me = 1840. 

7.1. The Electron Diffusion Scale 

The electrons have the same mean free path as the pro- 
tons, but they are faster than the protons by the square 
root of the mass ratio: 

C«,e«(^)'\, (57) 

where c^.e is the electron thermal speed and Cg is the 
sound speed, which is comparable to the proton thermal 
speed. Because of charge neutrality, the electrons have 
the same density, both perturbed and unperturbed, as the 
protons. Viscous damping caused by the electrons may 
be neglected; since the dynamic viscosity of the electrons 
is smaller than that of the protons by the square root of 
their mass ratio, electron viscous damping is always sub- 
dominant. 

As we now show, electrons are important for con- 
ducting heat on parallel lengthscales slightly larger than 
the proton mean free path. Electrons diffuse parallel 
to the magnetic field across a distance A|| in the time 
(Lmfp/cs.e)(A||/Lmfp)^. This is equal to the cascade time 
of an eddy with parallel lengthscale A||, i.e. it is equal to 
A||/w^, at the "electron diffusion scale", given by Ay ~ 
Led, where 

Led « ^L^fp « /3^/^(^)'^'Lmfp . (58) 

VA V me / 

In eddies with parallel lengths smaller than this, the elec- 
trons diffuse across the parallel lengths of the eddies, 
thereby conducting heat, and the electrons are isothermal; 
in eddies with parallel lengths larger than this, conduction 
is unimportant. 

We must also consider the effect of electron conduction 
on the protons. Electrons and proton temperatures ap- 
proach a common value on the timescale that there are 
mp/rrie collisions per particle. Since electrons are faster, 
the collision time is set by the electron speed, and the time 
for electron and proton temperatures to equilibrate is 

TOpL„ifp 1 1 /mp\i/2 

This time is smaller, by a factor /?, than the cascade time 
of eddies with parallel size Led. Thus in eddies of this par- 
allel size, and in those which are slightly smaller, protons 
are at the same temperature as the electrons. And, since 
the electrons are isothermal, so are the protons. 

As a result of the above considerations, when the par- 
allel cascade crosses the electron diffusion scale, the cas- 
cade becomes isothermal. This is similar to the crossing 
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of the cooling scale, discussed in §5.1, though backwards, 
and similar considerations apply. In particular, the Alfven 
mode is incompressible, and hence unaffected. The slow 
mode is nearly incompressible, and is only slightly affected: 
it suffers some damping to first order in 1/ p. The entropy 
mode, however, does not exist under isothermal condi- 
tions; as the electron diffusion scale is crossed, entropy 
waves are converted into slow waves, and the density fluc- 
tuations which had been associated with entropy waves 
now become associated with slow waves. 

7.2. The Electron-Proton Equilibration Scale 

Continuing to slightly smaller scales, the isothermal cas- 
cade reaches the "equilibration scale" , where the cascade 
time, A||/t;^, is comparable to the electron-proton equili- 
bration time (eq. [59]). At this scale. Am 



To evaluate the proton diffusion scale, we equate the 
time for protons to diffuse across an eddy of parallel 
lengthscale A|| , i.e. (Lmfp/cs)(A|| /Lmfp)^, with the cascade 
time, A||/i>A. This gives A|| ^ Lpd, where 



-lpd 



VA 



-^mfp 



/31/2L 



mfp 



(61) 



However, from an observational point of view, it is the 
transverse size of a damped eddy, , which would be ob- 
served. This is related to the parallel size through equa- 
tion (2). Alternatively, we use the Kolmogorov constant 
M|| = va)^^/v\^K\\ (see §2.5.1)."'^'^ Then, the transverse 
scale of a damped eddy, i.e. an eddy with parallel size 
A|| Lpd, is 

A±«M||^Lpd 



Leq, where 



VA 



1 /Wp\ 

i72 V^J 



1/2 



-'mfp 



(60) 



A I \i/3 



M,i/3L„ifp 



which we solve for the cut-off lengthscale: A_l 
where 



(62) 

t(-L) 



-(-L) 



M, 



3/2 



Lo 

1/2, 



VLout^ Vn/cm '*/ 



3/2 



cm 



On larger scales, the protons are thermally coupled to the 
electrons, and hence they are isothermal on slightly larger 
scales; on smaller scales, the protons are thermally inde- 
pendent of the electrons, and hence adiabatic. The tran- 
sition through the equilibration scale is nearly identical to 
the transition through the cooling scale (§5.1): the Alfven 
waves are unaffected, and they mix larger-scale isothermal 
slow waves into smaller-scale adiabatic slow waves and into 
smaller-scale entropy waves, although with some damping. 
Therefore, below the equilibration scale, the entropy mode 
reappears. 

The result of the calculations in both this subsection and 
the previous subsection is that, in crossing the electron 
diffusion scale and the equilibration scale, density fluctu- 
ations which were associated with the entropy mode on 
large scales are transferred from the entropy mode to the 
slow mode and then back to the entropy mode. Density 
fluctuations which were associated with the slow mode on 
large scales are unaffected.^ Therefore, for our purposes, 
these two lengthscalcs have little net effect on the density 
spectrum. Although there is some damping, the amount 
of damping is comparable to the amount at the cooling 
scale, and hence is not very significant. 

7.3. The Proton Diffusion Scale: Death of the Slow Mode 

and Entropy Mode 

Entropy waves and slow waves with parallel wavelengths 
smaller than the proton mean free path both damp on 
the timescale that protons stream across their parallel 
wavelengths (Barnes 1966). However, because we consider 
/3 > 1, there is a scale slightly larger than the proton mean 
free path at which these waves damp in the turbulent cas- 
cade. This is the proton diffusion scale, which is the scale 
at which protons can diffuse across an eddy within a cas- 
cade time. Below the proton diffusion scale, proton viscos- 
ity kills the slow waves and proton heat conduction kills 
the entropy waves. Alfven waves are unaffected by either 
of these effects. The density spectrum therefore cuts off 
below the proton diffusion scale. 

^ However, if the density fluctuations in the large-scale slow mode were less than the density fluctuations in the large-scale entropy mode, then, 
on small scales, the slow mode would be boosted so that its density fluctuations would be comparable to those of the entropy mode: the slow 
mode can steal approximately half of the entropy mode. 

Although we drop other factors of order unity, we keep the dependence on this one because of the importance of the proton diffusion scale. 



(63) 

, _3/ (64) 

^out^ \n/cm 

Below this scale, the density spectrum is cut off. For plau- 
sible values of /3, Lout, and n, this lengthscale is signifi- 
cantly smaller than L^fp. It is also larger than the proton 
gyroradius (eq. [66], below). 

Armstrong, Rickett, & Spangler (1995) summarize the 
observations of the density spectrum cut-off. There is con- 
siderable evidence that the cut-off scale is smaller than 
about 10^" cm along many directions. There is weaker ev- 
idence, from refractive scintillation, that along some lines 
of sight the cut-off scale is larger than around 10^ cm. Our 
theory might have implications for these observations. 

We conclude this section with two remarks. First, we 
re-emphasize the importance of the fact that the cascade 
is anisotropic. It is this fact which allows the density spec- 
trum to reach lengthscalcs which are significantly smaller 
than the proton mean free path. Second, we note an im- 
portant consequence of the paraUel cascades of the slow 
mode and of the passive scalar; these parallel cascades were 
explained in great detail when considering incompressible 
MHD turbulence (§2). Had there been no parallel cascade, 
then neither the slow mode nor the entropy mode would 
damp at the proton diffusion scale. The wavelength along 
the magnetic field would be effectively infinite, and so the 
ions could not diffuse across wavelengths. Rather, the slow 
mode and the entropy mode would be mixed down to the 
proton gyro scale, where the Alfvenic cascade is cut off. 

7.4. Density Spectrum Below the Proton Diffusion Scale 

On scales smaller than Lpj\ density fluctuations are 
wiped out. Since on these scales protons diffuse across the 
lengths of many eddies before the eddies cascade, the den- 
sity within neighbouring eddies is homogenized. Homog- 
enization occurs not only in the direction parallel to the 
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local mean magnetic field — as might have been expected 
since protons are tied to fieldlines — but also in the trans- 
verse direction. This is because eddies that are adjacent 
in the parallel direction incorporate substantially different 
fieldlines, so proton diffusion also wipes out density dif- 
ferences amongst eddies with transverse separations. The 
result is that the density spectrum on scales smaller than 
Lp^-* is determined by density fluctuations at Lp^^ ; equiv- 
alently, 

\i_ , for < lJ,-^^ . (65) 

8. THE END: THE PROTON GYRO SCALE 

The Alfven mode is undamped in a collisionless medium 
(Barnes 1979). Thus, the Alfven wave cascade survives be- 
low the proton diffusion scale, without the accompaniment 
of the slow waves and entropy waves. The Alfvenic cascade 
is damped at the scale of the proton gyroradius: 

Lp,gy, = V2cT^ = 2.5 X lO^f / cm . (66) 

eB \n/cm-^/ 

At this scale, Alfven waves arc converted into whistlers. 
The whistlers cascade to smaller scales, where they 
are damped by the collisionless effects of the electrons 
(Quataert 1998, viz. the curve in his Fig. lb that cor- 
responds to equal electron and proton temperatures). 

9. SUMMARY 

The primary goal of this paper has been to calculate the 
small-scale density spectrum in turbulent interstellar plas- 
mas. Our theory of compressible MHD turbulence is based 
upon the incompressible theory of Goldreich & Sridhar 
(1995), for which there is growing support from numerical 
simulations. We hypothesize that the compressible the- 
ory is similar to the incompressible theory of MHD tur- 
bulence, but with two main modifications: a compressible 
slow mode and an entropy mode which is passively ad- 
vectcd. While we believe that this hypothesis is plausible, 
it can, and should, be tested with numerical simulations 
of compressible MHD. 

Because of the multitude of special lengthscalcs encoun- 
tered when discussing the turbulent cascade in the inter- 
stellar medium, a recapitulation might be useful. In the 
following, we list the most important lengthscales, and 
summarize their significance to the turbulent cascade. 
The Outer Scale (Lout): This is the Icngthscale at which 
the turbulent motions are stirred. On slightly smaller 
scales, there is isothermal hydrodynamic turbulence for 
most plausible astrophysical sources. 
The MHD Scale (Lmhd): At this lengthscale there is 
a transition from hydrodynamic to MHD turbulence. 
Larger-scale hydrodynamic motions couple to smaller- 
scale Alfven waves and slow waves. 

The Cooling Scale (Lcooi): At this lengthscale there is a 
transition from isothermal to adiabatic turbulence. In 
high-/? turbulence, where small scale density fluctua- 
tions are due to the entropy mode, entropy — and hence 
density — fluctuations are suppressed by cooling. In low-/? 
turbulence, i.e. 1 ^ /3 < Lout/csicooi, density fluctuations 
due to the slow mode are important, and cooling has a 
negligible effect on the small-scale density spectrum. 
The Collisionless Scale of the Neutrals (Lh; Lne): Neu- 
trals decouple from ions across these scales. If the neutral 



fraction is not sufficiently small, then all three modes — 

Alfven, slow, and entropy - arc damped. 

The Collisionless Scale of the Ions (A|| = Lo^; Ay = Lcq; 

A II = Lpd X± = Lp^''): Across these lengthscales there 
is a gradual transition from fluid behaviour to collisionless 

plasma behaviour. The first two of these lengthscales, the 
electron diffusion scale and the equilibration scale, are set 
by the electrons. These lengthscales have only a small net 
effect on the cascade. The proton diffusion scale, however, 
is critically important for density fluctuations: both the 
slow mode and the entropy mode are cut off at this scale, 

and hence there are no density fluctuations below Lp'^^ 
The Proton Gyro Scale (Lp.gyr): Alfven waves are cut off 
at this scale. However, this scale has little importance for 
the density spectrum, because there are no density fluctu- 
ations below the proton diffusion scale. 

10. COMPARISON WITH HIGDON'S WORK 

In his 1984 paper, Higdon attributes small-scale density 

fluctuations to the passive mixing of the entropy mode. 
In his 1986 paper, Higdon notes that if the entropy mode 
varies along the magnetic field, it is damped in a collision- 
less medium. He then attributes small-scale density fluc- 
tuations to the passive mixing of tangential pressure bal- 
ances. Tangential pressure balances are structures which 
are parallel to the mean magnetic fleld. They are com- 
posed of both entropy waves and slow waves that have 
purely transverse wave vectors. 

Considering that Higdon's papers preceded even a the- 
ory of incompressible MHD tiu'bulence, they are a remark- 
able accomplishment. However, Higdon does not account 
for the parallel cascade. Entropy and slow waves with 
purely transverse wave vectors contain negligible power. 
They are cascaded along the magnetic field by Alfven 
waves. Consequently, they are damped when protons can 
diffuse across eddies in a cascade time, i.e. below the pro- 
ton diffusion scale. 

11. COMPRESSIBLE TURBULENCE WHEN /? < 1 

In this section only, we consider compressible turbulence 
in plasmas that have /3 < 1, e.g. in isothermal shocks and 
in the solar wind. 

Since Alfven waves are unaffected by the value of /3, 
and since nearly transverse slow waves arc only slightly 
affected (see Appendix), the dynamics of the cascade is 
nearly independent of /3. While the dispersion relation of 
the slow mode is changed from lv = VAkz to w = Cskz, the 
slow waves are still passively mixed by the Alfven waves. 
Therefore, the Alfven, slow, and entropy spectra in /? < 1 
MHD turbulence are the same as when /3 > 1. 

However, the damping of the slow and entropy waves is 
significantly changed. When /3 > 1, the proton thermal 
speed is faster than the Alfven speed. Therefore protons 
can stream across small eddies before they cascade. Con- 
versely, when /3 < 1, the time for protons to cross an eddy 
is always shorter than the cascade time. Therefore slow 
and entropy waves cannot be damped by protons which 
cross eddies, and the density spectrum extends to smaller 
scales. 

We have ignored the effects of the electrons. However, in 
the following we show that the electrons' behaviour may be 
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ignored for our purposes. Our discussion closely parallels 
that when /3 > 1 in §7, and uses similar notation. When 
(3 < 1, the equilibration scale (eq. [60]), is larger than 
the electron diffusion scale (eq. [58]). Thus the largest 
scale at which kinetic effects arc significant is the scale at 
which the cascade time is comparable to the time for elec- 
trons and protons to equilibrate their temperatures. Below 
this scale, electrons and protons are thermally decoupled. 
Nonetheless, this thermal decoupling has no effect on the 
Alfvcn, slow, or entropy waves on scales larger than the 
electron diffusion scale. At the electron diffusion scale, the 
cascade time is comparable to the electron diffusion time. 
Below this scale, the electrons are isothermal. However, 
since the protons are thermally decoupled from the elec- 
trons, the slow waves and entropy waves are cascaded to 
smaller scales. 

Continuing to smaller scales, the next scale of impor- 
tance is that at which the parallel size of an eddy is com- 
parable to the proton mean free path. Below this scale, the 
cascade is coUisionless. The Alfven waves are undamped 
by coUisionless effects. The entropy waves are undamped 
since protons cannot cross eddies within a cascade time. 



The slow waves are also undamped within a cascade time: 
although they damp within a waveperiod by Barnes damp- 
ing, their waveperiod is longer than the cascade time. Con- 
sequently, the Kolmogorov density spectrum extends down 
to the proton gyro scale. 

12. FUTURE WORK 

In a future paper, wo will examine; in detail the density 
spectrum in the solar wind. In another paper (paper II), 
we will relate the theory of compressible plasma turbu- 
lence developed here to observations primarily those of 
diffractive scintillation. We will demonstrate that the ob- 
served amplitude of small-scale density fluctuations is sur- 
prisingly large, especially along certain lines of sight, such 
as the one toward the Galactic Center. Then we will at- 
tempt to determine which astrophysical sources contribute 
the bulk of the density fluctuations and why they do so. 

We thank Jason Maron for analyzing some of his sim- 
ulations for us, for showing us how to use his numerical 
code, and for informative discussions. Research reported 
in this paper was supported by NSF grant 94-14232. 



APPENDIX 
WAVES IN COMPRESSIBLE MHD 



In this Appendix we derive the properties of the Alfven mode, the slow mode, and the entropy mode. The fast mode is 
not relevant to interstellar scintillation for reasons discussed in the body of the paper. The Fourier-transformed, linearized 
equations of ideal MHD, with d/dt ~iuj and V ik, read 

ton' — (fc • v)n , (mass conservation) (Al) 



uv ~ k 



c Vb^VA 



n 



bvAkz , (momentum conservation) (A2) 
ujb = zvA{k ■ v) — vvAkz , (Faraday's law) (A3) 

where n is number density, p is mass density, B is the background magnetic field intensity, and va = B/y^Anp is the 
Alfven speed. The z-axis is chosen to lie parallel to the background magnetic field. Some other variables are the perturbed 
number density, n', the fluid velocity, v, and the perturbed magnetic intensity, B' , or in velocity units, b = B'/y/Awp. 
The sound speed, c, is defined less conventionally here to be the square root of the ratio of perturbed pressure to perturbed 
mass density, c = (p'/p')^/^. For an isothermal gas. c = ct, where ct is the isothermal sound speed; = p/ p. For an 
adiabatic ideal gas, c = Cs, where Cs is the adiabatic sound speed; = jp/p = ^p/'ip for a monatomic plasma. More 
generally, however, the equation of energy conservation must be specified before c can be determined. While in general c 
depends on w, we ignore the w-dependence in this Appendix, and treat c as a constant. This treatment is valid for the 
purposes of this paper. 

We can express the equations of motion in terms of v as follows: 

(a;2 - klv\)v = [{k ■ v){c^ + v\) - A;^v^t;i]A; - v\kz{k ■ v)z . (A4) 
Alfven Mode 

The Alfven mode is incompressible; k ■ v = Q. Thus the term involving the soimd speed in the momentum equation 
vanishes, and the properties of the Alfven mode are independent of the equation of energy conservation. We obtain the 
dispersion relation 

uj = VA\kz\ , (A5) 
by forming the cross product of equation (A4) with k. The eigenfunction satisfies 

n' = k- v — z-v — , b — —sign{kz)v . (A6) 
Note that both v and b are perpendicular to z: the Alfven wave is polarized transverse to the unperturbed magnetic field. 
There is a change in the density spectrum at this scale which is of order unity; it is due to the change in the electrons' equation of state. 
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Slow Mode 

We summarize the properties of the slow mode to lowest order in fc^/fc <C 1, the limit appropriate to the MHD cascade. 
To obtain the dispersion relation wc assume, subject to verification, that the perturbation in total pressure i.e. thermal 
plus magnetic pressure — vanishes to second order in kz/k. Note that the perturbation in total pressure is proportional 
to the terms in square brackets in both equation (A2) and equation (A4). Prom the vanishing of this pressure term in 
equation (A4), we then have 

, 2n 



k ■ V = 



l + O 



{^)]kzv.. (A7) 



The z-component of equation (A4), with k ■ v given by equation (A7), then yields the dispersion relation: 



+ f A 



l + o(^^y]vA\k,\ . (A8) 



k 

Next, we solve for the eigenfunction to lowest order in kz/k: 

— ~ sign(A:^)^^^==L= , ~ - ;;2 , .,2 17^^ , t^y = , 6 ~ -sign(A;^)^==== , (A9) 

where the x-axis has been chosen to lie in the plane containing k and z. We now see that our assumption of negligible 

perturbed pressure is self-consistent to lowest order in k^/k. While wc have not used the x-component of equation (A4), 
both the left-hand side and the right-hand side of this equation are of the same order in kz/k. 

Note that both v and b are nearly parallel to z: the slow wave polarization is nearly aligned with the unperturbed 
magnetic field. Moreover, in the limit that c ^ va, the slow mode is nearly incompressible. 

Entropy Mode 

For adiabatic fluid motions, the linearized equation of energy conservation reads 

Lus' = , (AlO) 
where s is the entropy per particle. Thus, there exists a zero frequency, w = 0, mode with p' = v = b = 0, which has 

T- n " 7 5'' 
with T denoting temperature. Note that there is no corresponding mode in an isothermal fluid. 
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Fig. A1. — Wandering of Magnetic Fieldlines: a fieldline bundle of transverse size X± diverges after a parallel distance Ap where Ay is 
the parallel size of an eddy ( eq. [2]) as determined by critical balance. 
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Fig. A2. — Three- Dimensional Spectrum: contours of constant i'a^,A|| (sg- [3]), labelled by the value of v\j_^x^^/va- Along the X±-axis, 
■vXj_,o/'"A = (A_l/Lmhd)^/^; along the \\\-axis, vo,\^^/va = (Ah/Lmhd)^/^- 
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Table A1 
Summary of Lengthscales 

Symbol Description 



A_L lengthscale transverse to the mean magnetic field; 

it is comparable to the "observed" lengthscale. 

A|| lengthscale parallel to the mean magnetic field 

A|| parallel size of an eddy; it is a function of A_l (eq. [2]). 

Lout A_L = Lout at the outer scale 

Lmhd Ax = Lmhd at the MHD scale 

Lcooi Ax = Lcooi at the cooling scale 

Lh X± = Lh at the collisionless scale of hydrogen atoms 

Lne Ax — Lho at the collisionless scale of helium atoms 

Lmfp mean free path of protons and of electrons 

Led A II = Led at the electron diffusion scale 

Loq A II = Loq at the electron-proton equilibration scale 

Lpd A II = Lpd at the proton diffusion scale 

Lp^' A_L = Lp^' at the proton diffusion scale 

Lp,gyr A_L = Lp^gyr at the proton gyro scale 



